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Abstract 

Information to be stored and transported requires physical carriers. 
The quantum bit of information (qubit) can for instance be realised as 
the spin = degree of freedom of a massive particle like an electron or as 
the spin 1 polarisation of a massless photon. In this lecture, I first use 
irreducible representations of the rotation group to characterise the spin 
dynamics in a least redundant manner. Specifically, I describe the deco- 
herence dynamics of an arbitrary spin S coupled to a randomly fluctuating 
magnetic field in the Liouville space formalism. Secondly, I discuss the 
diffusive dynamics of the particle's position in space due to the presence 
of randomly placed impurities. Combining these two dynamics yields a 
coherent, unified picture of diffusive spin transport, as applicable to meso- 
scopic electronic devices or photons propagating in cold atomic clouds. 
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1 Introduction 



Classical information processing uses charge encoding where the bit values "0" 
and "1" are represented by a supplementary charge present or absent in a regis- 
ter. In the quantum limit, one eventually is led to consider a single elementary 
charge (electron or hole) in a quantum dot and hopes to realise quantum su- 
perpositions of "qubit" states |0) and |1) and entanglement between distinct 
qubits. But since they interact via long-range Coulomb forces, charge states 
suffer strongly from decoherence. Another discrete degree of freedom is spin. 
Spin i states, typically noted |±1) or |f , |), are the natural realisation of a qubit. 
But spin, and more generally information as such, needs a physical carrier. Can- 
didates here are electrons or holes (massive particle of spin s = i) and photons 
(massless, spin s = 1). A new promising field therefore is "spint(r)onics" : spin- 
based information transport and processing with electrons and photons. 

For these lectures, I chose to discuss the following model setting of diffusive 
spin transport (Fig. [T]): spin-polarised particles are injected from the left with 
probability pj (0) = 1 into a disordered sample and move diffusively towards 
the right, where a spin-sensitive detection reads out the final spin polarisation 
Pl(L) that we should calculate. 

The leitmotiv of this lecture was: Spin is a geometrical quantity, and one 
should be able to use irreducible representations of the rotation group in or- 
der to take advantage of symmetries. Since the participants of the school were 
not required to have been educated in group theory, I decided to give a rather 
complete introduction into representations of the rotation group with the hope 
that the sometimes intimidating jargon of group-theoretical arguments may be- 
come more familiar to everybody. Attending myself other lectures, I adopted a 
master-equation approach parallelling Klaus Hornberger's lecture [1] to which 
the present notes refer occasionally. During the lecture itself, only the first part 
on spin relaxation was delivered as presented in the following; in the second 
part on diffusion, I relied on arguments taken from diagrammatic perturbation 
theory that seemed to bewilder the audience more than anything else. The 
present notes remedy to this dissymmetry and treat also the diffusive part as 
momentum relaxation with a master equation. This parallel allows to combine 



p T (0) = l P T W=? 




L 

Figure 1 : Model setting of diffusive spin transport: spin-polarised particles are injected 
from the left with j>f (0) = 1 into a disordered sample and move diffusively towards the 
right, where a spin-sensitive detection reads out the spin-polarisation p-f (L) . 
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both dynamics rather economically into a single, coherent picture of diffusive 
spin dynamics; it is my hope that this conceptual unity is appreciated by my 
readers. 

These notes finish with a description of quantum corrections to diffusive 
spin transport including spin-flip effects. Other subjects covered during the 
school's lecture were of a more anecdotic type and, although hopefully enjoyed 
by the audience, did not seem to fit into the present format; readers interested 
in spintronics properly speaking are referred to the recent review by Zutic, 
Fabian, and Das Sarma [2J. I have to admit that I made no attempt to cover 
systematically the vast literature on the subject of irreversible spin dynamics 
and quantum transport, which would have been a hopeless task in any case; 
my apologies to many colleagues whose excellent contributions may not be duly 
cited in the following. 

2 Spin relaxation 

"Spin" is internal angular momentum [31 d] . This was recognised by Uhlenbeck 
and Goudsmit [5] following Pauli [6] who postulated the existence of a fourth 
quantum number to explain fine-structure features of atomic spectra. The clear- 
est experimental manifestation of quantised spin is arguably the Stern-Gerlach 
experiment [7 where silver atoms are deviated by an inhomogeneous magnetic 
field into two distinct spots on a detector screen. 

Dirac discovered that bispinors (vectors of four components, a spin | spinor 
and an anti-spinor) appear naturally when one looks for a Schrodinger-type 
wave equation in the relativistic framework of the four-dimensional Minkowski 
space. Wigner [8] showed that spin is one of the fundamental quantum numbers 
that permits to identify an elementary particle in the first place: it characterises 
the particle's properties under rotations in its proper rest frame. Therefore, to 
understand spin is to understand rotations. 

2.1 Spin — a primer on rotations 
2.1.1 Rotation group 

Physical objects are described by coordinates x — (x\, X2, £3) <E R 3 with respect 
to a reference frame in configuration space. Rotations (called "active" when 
the object is turned and "passive" when the reference frame is turned) are 
represented by 3x 3 matrices: x' = Rx. Proper rotations conserve the Euclidean 
scalar product x ■ y — ^\ Xiyi and the orientation of the frame. The rotation 
matrices are therefore members of SO (3), the set of orthogonal matrices RR t — 
R t R = I3 of unit determinant det R = +1. With the usual matrix multiplication 
as an internal composition law, theses matrices form a group, satisfying the 
group axioms: 

f . Internal composition: Vi?i, 2 £ SO(3) : i? 2 i = R2R1 £ SO(3); 

2. Existence of the identity: BE : RE = ER = R MR £ SO(3) with E = 1 3 ; 
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3. Existence of the inverse: VRBR' 1 : RRr 1 = R~ X R = E. 



The group is non-Abelian because the matrices do not commute: R2R1 ^ R\R%. 
An exception are rotations of the plane around one and the same axis, forming 
the Abelian group SO (2). 

A possible parametrisation of a rotation is the polar description (n, 6) =: 6 
with h the unit vector along the rotation axis and 9 £ [0, 7r] the rotation angle. 
The following two rotations of configuration space are identical: 

i?(n, = tt) = R(-n, 9 = tt). (1) 

These opposite points must be identified such that there are closed parame- 
ter curves that cannot be contracted into a single point, see figure [21 This 
means that SO (3) is a doubly connected manifold. Instead of studying this 
projective group, one may also turn to its universal covering group SU(2), the 
group of all unitary 2x2 matrices over (D with unit determinant. SU(2) is sim- 
ply connected, and there is a homomorphism (a mapping preserving the group 
structure) linking every element U £ SU(2) to the rotation R £ SO(3): the 
rotation x' = R(6)x is described by 

x' ■ a = U(0)x ■ <tU{0)\ (2) 

Here, er = (a\, o~2, 03) is a vector whose entries Ci are the Pauli matrices such 
that 

x-*=( Xl ~' 1X A (3) 

\x ± + ix 2 -x 3 ) 

and the unitary rotation matrix acting from the left and from the right is given 

by 

9 9 , • cr . 

U(9) = 1 2 cos — — m ■ cr sin - = exp(— 1 — - — ). (4) 



Figure 2: The parameter space of the rotation group SO(3), a filled sphere of radius 
7r. Two identical rotations [eq. {!}/ can be connected by a closed curve that is not 
reducible to a single point. 
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Note that this mapping is two-to-one because U and — U yield the same rotation, 
and that it takes a rotation by an angle of Air to recover the identity transfor- 
mation: U(9 = 2ir) = — 1 2 , but U(6 = Air) = +I2. This is not a mysterious 
quantum property as is sometimes stated, but reflects the double-connectedness 
of SO(3)-rotations in our everyday reference framed Dirac's construction of a 
solid body connected by strings to a reference frame is supposed to convey an 
"experimental" idea of this property [4]. 

2.1.2 Representations 

A group G can act in many different disguises that share the same abstract group 
structure, as defined by the multiplication law or group table. These different 
appearances are called (linear) representations. Mathematically, they are map- 
pings T>i : G — > GL(Vi) from the group G to the general (linear) group of regular 
transformations D : Vj — ► Vi of a vector space Vi into itselfd Importantly, this 
mapping must be a homomorphism which means that the representation has the 
same group structure as G. Notably, for all elements 31 , (72 £ G (with the prod- 
uct <?2ffi S G) the representing transformations verify D{g 2 g\) = D{g2)D{g\). 

AmiVi is called the dimension of the representation T>i. Finite-dimensional 
linear representations are given in terms of quadratic matrices of size dim Vi . A 
representation is called reducible if there is a basis of V = Vi ® V2 such that all 
transformations D 6 Z> arc written 

°={ D o i)- («) 

In other words: the transformations D\ : V\ — > V\ define already a representa- 
tion T>\ of its own. If the matrix is block-diagonal (* = 0), the representation is 
completely reducible. If it is not reducible, it is called irreducible^ meaning that 
one has achieved to work in the smallest possible subspace. 

The rotation group SU(2) is a manifold that depends on a continuous set of 
parameters 6 = (01,6*2,6*3) with respect to which it is infinitely differentiable. 
This structure is called a Lie group. The group multiplication law is completely 
determined by the commutation relation of its generators: 

[Jj, Jk] ■= JjJk — JkJj = itiejkiJi- (6) 

Here as in the following, the sum over repeated indices is understood. These 
generators are said to form a (representation of) the Lie algebra su(2). In the 
so-called natural representation of SU(2) by itself [eq. (0])], the generators are 
J = Other representations will feature different generators, but all group 
representations share the same commutation relation! Finite rotations are gener- 
ated by exponentiation: U{8) = exp(— iO ■ J /ft). Topologically speaking, SU(2) 
is compact. According to a general theorem, all representions of a compact Lie 
group are completely reducible to finite-dimensional irreducible representations. 

1 Cartan developed half-integer spin representations as early as 1913 within the theory of 
projective groups [9]. 

2 T> stands for the German word "Darstellung" ; in anglo-saxon literature, often the symbol 
r is used. 
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2.1.3 Functional representation of SO(3) 



As an example for an infinite-dimensional representation of the rotation group 
S0(3) already useful in classical physics, consider the transformation of func- 
tions / describing the position of an object on the unit sphere. Saying that the 
object is rotated to r' = Rr implies that the function is transformed by 



The action of D in the (infinite-dimensional) functional space can be written 
D = exp{— i9-T} where T = — ir x V is a differential operator, as can be verified 
by considering an infinitesimal rotation around 9 = 8n: to first order in 8, we 



have /(iTV) = f(r - 6 x r) = f(r) - {6 X r) ■ V/(r) = [1 - 6 ■ (r x V)]/(r). 



The finite-dimensional irreducible representations to which this infinite-di- 
mensional one can be reduced are obtained by decomposing / into surface har- 
monics Yhm\ each subspace L = 0, 1, . . . then admits an irreducible represen- 
tation of dimension 2L + 1. Surface harmonics are a concept arising already 
with the multipole expansion of charge distributions in classical electrodynam- 
ics. But to cite Hermann Weyl [10]: "This reveals the true significance of surface 
harmonics; they are characterised by the fundamental symmetry properties here 
developed, and the solution of the potential equation in polar co-ordinates is 
merely an accidental approach to their theory." 

2.1.4 What is quantum about spin? 

Spin, and especially half-integer spin, is much less mysterious than sometimes 
suggested by standard textbook wisdom. Half-integer spin does not require the 
kinematic framework of special relativity, it arises already in Galilean relativity. 
Also, as mentioned at the end of l2.1.1[ reference frames for solid bodies already 
introduce half-integer spin. However, there are genuine quantum features to be 
aware of: 

1. In classical mechanics, Noether's theorem assures that if the Lagrangian 
function is invariant under infinitesimal rotations, then the orbital angular 
momentum L is a conserved quantity. However, L a priori has nothing 
to do with the generator of rotations T introduced in 12. 1 .31 In quantum 
mechanics, thanks to the appearance of h, the generators can be identified 
with observables J — TiT with dimension of angular momentum. If the 
Hamiltonian is invariant under all rotations U(0), their generators itself 
are conserved quantities, and Noether's theorem takes a very simple form: 



Separating the orbital part L of angular momentum from the total angular 
momentum J = L ® Is + Ix, ® S =: L + S, one identifies the rest-frame 
angular momentum or spin S obeying the same fundamental commutation 
relations (JHJ). In the remainder, we will only have to consider the spin part 
and let L = 0. 



(Df)(r) = /(iTV). 



(7) 



H' = UHU^ [H, U] = & [H, J] = j = 0. 



(8) 
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2. The observables S generate irreducible representations of dimension 
d s = 2s+l with s = 0, o> • • • and discrete magnetic quantum numbers m — 
— s, — s + 1, . . . , s. Pure states are noted |sm). The Casimir operator S 2 
specifies the irreducible representation, S 2 \sm) = h 2 s(s + l)\sm), whereas 
the magnetic quantum number gives the projection of the spin onto the 
quantisation axis (usually called the z-axis): S z \sm) = Tim\sm). 

3. States can be classified regarding their transformation properties, but here 
with a more general importance than in classical mechanics due to the 
superposition principle. An atomic s-orbital, for instance, may be seen as 
a "coherent superposition of all possible Kepler orbits" and is invariant 
under all rotations. In section I2.3.3) we will introduce the irreducible 
components of mixed quantum states, also known as state multipoles. 

2.2 Master equation approach to spin relaxation 

As a simple model for spin dynamics, we shall study the Hamiltonian 

H = -fj,S B . (9) 

It describes the coupling of the magnetic moment fi = pS to a magnetic field 
B. For electrons, p — — gpB in terms of the Bohr magneton /zb = \e\h/(2m e c) 
and the gyromagnetic ratio g — 2.003... in vacuum. The spin operator has 
been chosen dimensionless such that its action in the irreducible representation 
2?( s ) will be S z \sm) = m\sm) and S 2 \sm) = s(s + l)\sm) for the remainder 
of the lecture. The density matrix or statistical operator p of a spin S is a 
positive linear operator of trace unity on the state Hilbert space Ti s — <D rfs of 
dimension d s := dim7Y s = 2s + 1 that determines the expectation values of 
abitrary observables O as (O) = tr{pO}. 

2.2.1 Unitary spin dynamics 

According to one of the fundamental axioms of quantum theory, any closed 
quantum system evolves unitarily according to the Liouville-von Neumann equa- 
tion 

ihd tP = [H,p]. (10) 

This equation of motion is formally solved as pfo) — U(t2,ti)p(ti)U(t2,tiy by 
applying the time evolution operator for a time-dependent Hamiltonian, 

C/(t 2 ,t 1 ) = Texp|-^ 3 J H-(t')dt / | (11) 

where T[fr(ii)fl"(t 2 ) . . . H{t n )\ = H^Hfa) . . . H{t k ) for U > tj > ■ ■ ■ > t k is 
Dyson's time-ordering operation. 
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2.2.2 Non-unitary spin dynamics: a classical model derivation 

Phenomena like "relaxation, damping, dephasing, decoherence,..." have in com- 
mon irreversible dynamics with an "arrow of time" [11) due to the irrevocable 
loss of energy and/or information into inobservable degrees of freedom which 
are usually called "bath" or "environment" . 

As an introductory model, we consider the Hamiltonian H = — [iS ■ B(t) 
with a randomly fluctuating magnetic field B(t). Predictions about the spin 
will involve an average over the field fluctuations which we describe as a clas- 
sical stochastic process [12] . This approach is typical for the physics of nuclear 
magnetic resonance; a regular driving magnetic field can of course be included 
in the treatment, but here we concentrate on the effect of random fluctuations. 

To obtain the effective dynamics, we develop the time-propagated density 
matrix p(t+ At) = U(t + At, t)p(t)U(t-\-lS.t 1 ty to second order in the interaction 
Hamiltonian, 

t+At t+At 

p(t + At)=p(t)-± J [H(h),p(t)]dh + J J H(t 1 )p{t)H(t 2 )&t 1 &t 2 
t t 

t+At 

~~^JJ T ( fr (*i)^(^)p(t)+p(t)i?(ti)ff(t2))dtidt 2 , 



t 



(12) 

and perform an average over all possible realisations of the fluctuating mag- 
netic field B(t). As a stochastic process [12], it is completely specified by its 
correlation functions 



Ciii2...i n (tl,t2, ■■■,tn) '■= B il (ti)B i2 (t 2 ) ■ ■ .B iti (t n ), (13) 

the ovcrline indicating an ensemble average over the field distribution. This 
distribution is taken to be centered Gaussian such that all correlation func- 
tions factorise into products of pair correlations. Therefore, only the first two 
moments need to be specified: Bi{t) = (zero mean) and 

Bi(ti)Bj(t2) =: B*Cij(ti,t2). (14) 

We assume a stationary process that depends only on the time difference ti — t 2 
and has a very short internal correlation time r c such that (Hj{t) = CijT c S(t). 
This last assumption of "white noise" (the power spectrum Cij{u) = est. contains 
all frequencies with equal weight) is valid if the noise correlation time r c is much 
shorter than the relevant time scale of the spin dynamics that is still to be 
determined. Lastly, we assume that the fluctuations are isotropic, Cij — ^Sij. 

Now we average the time-propagated density matrix (|12p over the field fluc- 
tuations and use the Born assumption that there are no further correlations 
between the fields appearing explicitly and the average density matrix. Then, 
the ensemble average applies to the fields only. The term linear in H disappears 
because Bi = 0. In the second order terms, one of the time integrations is con- 
tracted by the 8{t\ — ^-distribution of the correlation function; the remaining 
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integrand is time independent such that the integral gives just a factor At. The 
average time-evolved density matrix then becomes 

p(t + At) - (1 - 7.Ai)p(t) + 7s At &p(*)& + 0((jsAt) 2 ) (15) 

i 

where S, := Si/y/s(s + l) is the "normalised spin operator" with J2 i S 2 = 1. 
The spin relaxation rate 7 S := s(s + \)uj 2 t c is given in terms of the squared 
effective Larmor frequency lo 2 — p 2 B 2 /(3h 2 ). The relevant time scale of evolu- 
tion turns out to be t s :— 1/7,5. The effective time-evolution (|15p then is valid 
for a small time step At <C t s such that indeed 7,5 At <C 1. 

Exercise 1 - Non-unitary spin dynamics: quantum derivation 

Consider the time-independent model Hamiltonian H = —TiJS ■ t where our spin 
S is coupled to a freely orientable magnetic impurity, here modeled as a spin | with 
Pauli matrices t. The effective spin dynamics of S is described by its reduced density 
matrix p(t) = tr r {ps-r(t)} obtained by tracing out the uncontrolled impurity spin. 
Develop the time-evolved complete density matrix psv(t) to second order in J as in 
the previous section 12.2.21 and take the trace over the impurity with initial statistical 
mixture p T (0) = 5I2 (it is helpful to use the identities tr T {ri} = and tr T {rjTj} = 
2Sij). Show that one finds exactly the evolution (|15(l with a relaxation rate 7 S = 
s(s + l)J 2 At up to higher-order terms that become negligible in the formal limit 
At — > together with J — > 00 taken at constant 7 S (f| 



2.2.3 The quantum channel and its operator sum representation 

In the language of quantum information, the time evolution (|15|) up to order At 
defines a quantum channel p — p(t) 1— ► p' = p(t + At) and is here given in the 
so-called operator sum representation [131 8.2.3.] 

3 

p' = Y< w *p w i ( 16 ) 

with 

W := Vl- 7s Atl, Wi:=V%MSi, i= 1,2,3. (17) 

It is easy to verify that X^=o WiW} = 1 which guarantees the trace conservation 
trp' = trp. Kraus has proved that a channel of this form assures that the final 
density matrix is again completely positive. Therefore, it is also known as the 
Kraus representation and the Wi are commonly referred to as Kraus operators, 
cf. sec. 3.1 in [1]. In contrast to the unitary evolution p' = UpU^ of (|10|) . the 

3 But attention: at finite coupling J, the dynamics of our spin shows recurrence on a time 
scale given by the so-called Poincare time t TCC oc l/J. One could obtain a truly irreversible 
dynamics only by supposing that the single impurity spin is reset rapidly enough in order to 
dispose of the coherence. Alternatively, one may imagine the setting treated in section [4. 1.21 
our spin is moving and encounters different impurity spins such that in the thermodynamic 
limit the Poincare recurrence time goes to infinity and true irreversibility sets in. 
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appearance of several independent terms in the sum (|16[) signals non-unitary 
dynamics. 

For a spin i with Si = Cj/ V3, this quantum channel is the qubit depolarising 
channel |13[ 8.3.4.]: 

/ ) ' = ( 1 -Pi)/'+yX! (7 ''* T " Pi = 7sAi. (18) 

With equal probability pi/3, the qubit is affected by the action of one of the 
Pauli matrices <7j, and with probability 1 — p±, it remains untouched. Since for 
spin | one may write ^ Oipcn — 23.2 — p, the channel also takes the suggestive 
form 

p ' = (1- P2 ) p+ Plt 2 . (19) 

This means that with probability p 2 = ^P\/3 (remember pi <C 1 such that 
also P2 <C 1), the density matrix is taken to a complete statistical mixture and 
remains identical with probability (1 — p 2 ). 

By convention, the "depolarising channel" for higher spin S > 1 (i.e. d s x d s 
density matrices with d s — 2s + 1 > 3) is still defined via the relation (fT9| with 
the statistical mixture ;j-ld 3 on the right-hand side. This channel is also called 
"SU(n) channel" with n = d s because the corresponding Kraus operators are the 
n 2 — 1 generators of the Lie algebra su(n) [14]. Note that our physical model of 
an arbitrary spin S coupled to a fluctuating magnetic field does not lead to this 
specific Lie algebra channel: obviously, the operator-sum representation (I16|) 
contains only the 3 generators of su(2), albeit in a representation of dimension 
d s = 2.s + 1. We will see in the following that this makes the spin dynamics 
richer and its description more involved. Group-theoretical methods will be 
introduced that are well adapted to cope with this complexity. 

2.2.4 The Liouvillian 

The linear operators on the state Hilbert space Ti s are themselves elements 
of a linear vector space (we can add operators and multiply them by complex 
numbers). This vector space is called Liouville space L{TL S ) and is spanned, for 
example, by the basis of dyadics induced by basis vectors \n) of TL S : 

|m)(n| =: |t7ot), n, m = 1, . . . , d s . (20) 

The Liouville- von Neumann equation of motion idtp = Cp for a closed quantum 
system defines the Liouvillian 

t = \[ H i-] (21) 
whose matrix elements in the dyadic basis are 

( 7777?, | £ 1 772 71 ) — £"mn,m'n' — T" {^H nirn i 6 nn > H n ' n Smm'^) • (^^) 
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Exercise 2 - Liouvillian eigenvalues 

Show that the eigenvalues of the Liouvillian in the basis {|mn)} induced by the 
energy basis H\m) — e m \m) are the possible transition frequencies uj mn — (e m — 
e n )/h. These are experimentally accessible quantities, in contrast to the absolute 
energy eigenvalues of the Hamiltonian H. 

As an operator between operators, the Liouvillian C : L(Tt s ) — > L(H S ) is 
called a superoperator jT^l HE] El HE] • The superoperator formalism is a con- 
venient starting point for projection operator techniques in statistical dynamics 
[16] , effective dynamics of open quantum systems [19] and time-dependent per- 
turbation theory [20] . 

2.2.5 The Lindbladian 

The effective dynamics of our spin coupled to a randomly fluctuating field 
can also be formulated in terms of a superoperator. Taking the formal limit 
liniAt-^o p ( t+z ^j £ffl — : dtp(t) in (|15[) leads to the master equation 



a linear equation of motion for the effective density operator, whose effective 
generator of time evolution is the Lindbladian [18] 



= - y & /"(*)]] = |E (2&p(*)& - SiS iP (t) - P (t)SiSi) . 



This is a pure relaxation superoperator in the Lindblad form, which assures 
the complete positivity of the time-evolved density matrix. The symbolic limit 
At — > really means At'jg <C 1 but is still assumed to be "coarse-grained" 
compared to the field correlations, At 3> t c . In this limit, the master equation 
(|23[) is a linear equation for the density matrix p(t) local in time t and thus 
describes Markovian dynamics without any memory effects. 

Exercise 3 - Lindbladian matrix elements 

Show that the superoperator matrix elements of the Lindbladian (|24[) are given by 



and verify the trace-preserving property C mm ,m'n' = from this expression. 



2.3 Irreducible scalar spin relaxation rates 

Formally, the solution of the master equation (|23[) is very simple: 



d tP (t) =£p(t), 



(23) 



(24) 




(25) 




(26) 
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The dynamics induced by the Lindbladian is called a "quantum dynamical semi- 
group" (Sec. 3.1 in 1 ]). Indeed, the time evolution superoperator satisfies 
cxp[C(t 2 + ti)] = exp[£t 2 ] exp[££i] and exp[£0] = 1 which indicates a group 
structure. "Semigroup" means that the inverse to each element does not need 
to exist, and indeed here it does not since the non-unitary dynamics obtained 
by tracing out the environment has induced an arrow of time. The Kraus repre- 
sentation (fT6|) and the Lindblad form ([24)) guarantee the complete positivity of 
the final density matrix if the initial one is completely positive, but the inverse 
is not true: the Lindbladian is not invertible [HI 3.4.1.]. We will indeed see in 
section [2.3.51 below that C has one vanishing eigenvalue Ao = 0. 

In Liouville space, the master equation is a matrix equation dt\p) = £\p) 
of dimension d 2 x d 2 . As always when dealing with matrix equations, we have 
to diagonalise the Lindbladian in Liouville space in order to be able to use 
the formal solution (f26|) . For spin \ and d s — 2, diagonalizing a 4 x 4 matix 
is elementary and can be done by hand, but already for spin 1 with d 2 s = 9 
this becomes cumbersome. For higher spin, one definitely needs to resort to an 
efficient strategy to find the eigenstructure of the Lindbladian. In the following, 
we show how to take maximal advantage of rotational symmetries by using 
group-theoretical methods that lead to a very simple and physically transparent 
description for the effective spin dynamics. 



2.3.1 Scalar relaxation process: what results we should expect 



The Lindbladian was obtained by an isotropic average and is thus a scalar 
object, i.e., invariant under rotations. A rather high-brow way of expressing this 
simple property is to say "it transforms under the trivial representation X>(°)" . 
We may anticipate that the statistical operator can be decomposed into parts 
that transform under the irreducible representations T>^ K ^ of the rotation group, 
p = Y, K P (K) - The Lindbladian as a scalar object can only connect subspaces of 
equal rank K. Furthermore, inside each subspace it cannot distinguish between 
different orientations. Thus, in an adapted basis of Liouville space, it can be 
written as a purely diagonal matrix 

/Ao . . . \ 
Ai ... 



C = 



V 



(27) 



A<j|/ 



How many different eigenvalues may we expect? The total number is the di- 
mension of the Liouville space, d 2 s . Each subspace of rank K will have dimension 
dx = 2if + 1. Therefore, we will have to find only 2s + 1 different eigenvalues 
A x , each of which has degeneracy 2K + 1: 



particle 


s 


d s 


dl 


eigenvalues 


degeneracy 


electron 


l 

2 


2 


4 


Ao, Ai 


1, 3 


photon 


1 


3 


9 


Ao, Ai, A2 


1, 3, 5 
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Can some of the eigenvalues be identical which would imply an even larger 
degeneracy? A trivial example for this would be any operator proportional to the 
identity. But we will see below that this is not the case for the Lindbladian (f!M|) : 
the eigenvalues pertaining to different subspaces K K' are indeed different, 
Ak 7^ Aj<-'- We are therefore sure to have reduced the problem to the simplest 
possible formulation. 

2.3.2 Irreducible tensor operators 

Before we can define irreducible superoperators, we first had better understand 
the simpler concept of ordinary irreducible operators. An irreducible tensor 
operator by definition is a set of 2K + 1 components Tq, Q = —K, —K + 
1, . . . ,K, that transform under irreducible representations of rank K (i.e., whose 
transformation does not mix different K): 

{T§)' = UT§rf= J2 D™T$. (28) 
Q=-K 

Equivalently, one specifies the infinitesimal rotation properties by requiring 

[J±,T$\ = h^K(K + l)-Q{Q±l)T ( Q % (29) 
[J ,T$] = hQT*. (30) 

Here, the angular momentum raising and lowering operators are J± = J x ± iJ v . 
The simplest examples for irreducible tensor operators we need to know for the 
following are 

- K — or scalar operator Tq, a single operator that commutes with all com- 
ponents of the total angular momentum J. For instance S 2 , the Casimir 
operator indexing the irreducible representations T>^ . 

- K = 1 or vector operator with Euclidean components A = (Ai, A\, A3) 
satisfying 

[Jj,A k ]=me jk iAi, (31) 
and spherical components ^4o = ^.3 and A±\ = ^f-k=(Ai ± 1A2). 



Exercise 4 - Irreducible or not? 

2 

Is the Hamiltonian H = 2— — /j,S • B of a free massive particle coupled to an 
external magnetic field B via the magnetic moment ^1 = /j,S of its spin S a scalar? 
An irreducible tensor operator? 
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2.3.3 State multipoles 

In the usual Hilbert space basis, the statistical operator reads 

P= ^ Pmm'\sm){sm!\. (32) 

m,m' 

Here, the ket \sm) transforms under the irreducible representation T)( s > , whereas 
the bra (sm'|, as its complex conjugate, transforms under (V^)* , the contra- 
gredient representation. The ket-bra \sm)(sm / \ transforms under the direct 
product X>W (T>( s >)* which is reducible. One introduces therefore an ensem- 
ble of elements that do transform under the irreducible representation T>^ K \ 

Tq :=T${s,s) := ^ (-) s -" l (W -m\KQ) \sm!) (sm\, (33) 

m,m f 

with K = 0, 1, . . . , 2s and Q = -K, -K+1,...,K. The Clebsch-Gordan coeffi- 
cients {siS2mim2\KQ) are the coefficients of the unitary basis change from the 
direct product H S1 ®'H S1 towards the Hilbert subspace Ti. K of spin K that should 
be familiar from the addition of two spins. The CG coefficients are non-zero 
only if two selection rules are satisfied: (i) the two magnetic quantum numbers 
on the left add up to the one on the right, mi + ra-i — Q\ (ii) the angular 
momentum on the right satisfies the triangle inquality \s\ — S2I < K < s% + S2- 

In our case, we do not couple two spins, but a spin and its complex conjugate. 
Since there is a contragredient representation in the game, the CG coefficients 
in (|33p feature a characteristic minus sign in front of — m. The triangle selection 
rule implies for us < K < 2s. 

The irreducible tensor operators Tq =: \KQ) form a basis of Liouville space 
that is properly orthonormal with respect to the trace scalar product of matrices: 

(KQ\K'Q') := tr{(T#)tT#'} = <W<W- (34) 

The Hermitian conjugate is {Tq)^ = (—)QT^q. Any linear operator A can be 
decomposed in this basis, 

A = Y,A KQ T Q X with A KQ ;={KQ\A) = tv{(T^A}. (35) 

KQ 

The irreducible components 

PKQ = tr{(T*)tp} = ((T*)t) (36) 

of the density matrix are called state multipoles or statistical tensors, and have 
been introduced already in the 1950's by Fano and Racah [2"Tl I22j . 

Exercise 5 - Irreducible tensor operators 

(0) Show that all Tq except T ° have zero trace and calculate the state monopole 
moment poo (use \/2s + l{ssm'—m\00) = (— ) s ~ m 5 mrn i). 
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(1) Show that the irreducible vector operator is proportional to the spin operator, 
Tq — cJ 1 ^ 2 Sq. Hint: consider the action of (y / c7(ssm — raj 10} = (— ) s ~ m m), argue 
with rotational invariance and fix the proportionality constant c s = s(s + l)d s /3 by 
computing tr[(T 1 ) 2 ]. 



2.3.4 Irreducible spin superoperators 

By inserting the decomposition (|35[) on the left and right side of an arbitrary 
superoperator acting like CA, one obtains 

LA= \K'Q'){K'Q'\C\KQ)(KQ\A) (37) 

KQK'Q' 

such that with the notation Ck'Q',kq ■— (K'Q'\C\K Q) the superoperator reads 
C= Y, £k>q>,kq\K'Q')(KQ\ (38) 

KQK'Q' 

where the Liouville-space dyadics on the right hand side transform under T>( K 

, in an analogous manner to f|32|) . Following the same strategy as previ- 
ously, we recouple the elements again using the appropriate CG-coefficients to 
get irreducible superoperators [23j 

T&(K,K') := Y {-) K - Q (K'KQ'-Q\LM) \K'Q'){KQ\ . (39) 
Q,Q' 

Finally, any superoperator in completely decomposed form reads 

£ = EE C ™( K > K')TUK, K>) (40) 

LM K,K' 

with coefficients Clm(K, K') = Y,q, Q '(-) K ~ Q (K'KQ'-Q\LM)£ k , q ,, kq . This 
decomposition is completely general and applies to arbitrary superoperators. It 
is only worth the effort, however, if the superoperator has rotational symme- 
tries. The greatest gain in computational speed and conceptual clarity is ob- 
tained if the superoperator is a scalar such that its only nonzero component is 
L = 0,M = 0. In that case, which applies to our Lindbladian ([24]) . one finds 
by virtue of the triangle rule that K — K'\ as promised, C indeed connects 
subspaces of equal rank. We can choose the decomposition 

C=Y ^kT (K) (41) 

K=0 

where the T^ K ^ are orthogonal projectors onto the subspaces L(Ti s )^- K ' > of irre- 
ducible tensor operators of rank K: 

TW = V2K + 1T °(K, K) = Y \KQ)(KQ\ . (42) 

Q 
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They are indeed orthogonal, T^T^ = 8 KKI T^ K \ by virtue of I® and 
resolve the identity in Liouville space, ^2 K T^ K ' = 1, by virtue of a completeness 
relation of CG-coefhcients. 

2.3.5 Scalar relaxation rates 

Once the invariant subspaces are known, the eigenvalues are obtained by pro- 
jecting the superoperator onto an arbitrary basis element: 

\ K = (KQ\£\KQ) = tr {(T^CT«) . (43) 

For the spin relaxation Lindbladian (|24[) , one must calculate commutators of 
the form [S», [S^Tq]]. In order to use the defining commutation relations (|29j) 
and (|30|) for irreducible tensor operators, one writes the scalar product of spin 
operators in terms of the spherical components Sq — S z , S±i — ~F^(Sx ± iS y ), 

^2SiSi= Yl ( S p) fS P= E (-TS P S- P (44) 

i p=0,±l p=0,±l 

such that 

XK = -2£TT) S H P tr{{T^[S- P ,[S P ,T§]]} (45) 

It is now a simple exercise to show with the help of (|29|) and (f30]) (by paying at- 
tention to the supplementary factor \f2Ti in the definition of the raising/lowering 
components J± — \/2hS±i) that the double commutator gives back the tensor 
operator itself, [S- p , [S p ,Tq]] — c(p 7 Q, K)Tq . Summing all three terms gives 
remarkably simple eigenvalues, 

K(K+1) , N 

\k = ~1 s \ K = 0,1,... ,2s. 46 

2s(s + 1) 

These real, negative eigenvalues describe a pure relaxation process as expected 
from the definition of the Lindbladian (|24|) . They are scalar objects, also known 
as rotational invariants, and can be expressed in terms of 6j-coefficients that are 
constructed out of the irreducible representations T>( K > for the state multipole, 
X>( s ) for the spin S itself and T>^ for its coupling to the vector field B [A, 
123) . This type of consideration is of considerable importance in many different 
fields of physics involving angular momentum or spin; for example, relaxation 
coefficients very similar to (|46[) characterise spatial correlations in certain ground 
states of quantum spin chains [24j . 

2.3.6 Isotropic spin relaxation 

The master equation dtp(t) = Cp{i) separates into uncoupled equations for each 
invariant subspace: 

d tP (t) = Y, dtPK Q (i)T« = PkqWT* = J2 PKQ(t)\ K T« (47) 

KQ KQ KQ 
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The resulting differential equation d t pKQ{t) = ^KpKQ(t) for the state multi- 
poles is easily solved to yield a simple exponential decay 

PK Q (t) = e-^Wg(O) (48) 

with relaxation rates = | Ajf| . This is a particular example for a state multi- 
pole relaxation as described by Blum in chapter 8 of [ST]. The first two values 
deserve a special discussion. The scalar mode relaxation rate 70 = assures 
the trace preservation of p(t). At the same time, this vanishing eigenvalue is 
responsible for the fact that the Lindbladian is not invertible. The vector mode 
relaxation rate 71 = j s /[s(s + 1)] =: 1/n describes the relaxation of the orien- 
tation or average spin vector: 

(S(t)) = tr{p(t)S] = e-*M (5(0)) (49) 

since tr{Tg Sq>} projects onto K = 1 (remember exercise 5.(1) and the orthog- 
onality relation (f54"| ). 

For a qubit spin ^ this is all that needs to be calculated since any 2x2 
density matrix can be parametrised as p%{t) = |la + (S(t)) ■ a and thus 

p 2 (*) = e-*/>(0) + (l-e-^)il 2 . (50) 

This isotropic spin ^ relaxation therefore is for all times given by the depolarising 
channel (fl"9)) with p%{t) = 1 — e~'/ Tl . The Kraus operators for the operator sum 
representation (fit))) valid for all times follow by using p x (t) — jp 2 (t) : 



Vl + 3c-'/ T i VI - e _t / Tl 
Wo(t) = 1 2) Wi(t) = en, z = 1,2, 3. (51) 

Naturally, by developing these expressions to first order in tjr\ = ^"f s t one 
finds the Kraus operators for an infinitesimal time step derived in section !^. 2. 31 
In general, it is easy to show by derivation that to each quantum dynami- 
cal semigroup described by an exponential superoperator exp[£<] corresponds 
a Lindblad-type master equation [35], cf. Sec. 3. 2. in pQ. However, as always, 
the inverse operation of integrating the infinitesimal time evolution to finite 
times is much harder. Deriving a set of Kraus operators for an arbitrary quan- 
tum channel in general requires the complete diagonalisation of the microscopic 
Hamiltonian. Luckily, for spin ^ everything is so simple that the complete 
calculation is possible. 

Naturally, one thus wonders whether the full-fledged angular momentum 
formalism is necessary at all to describe isotropic spin relaxation. A simple 
calculation shows that in the "depolarising channel" defined by (fT9|h all non- 
scalar eigenvalues Ai, A2, • • • are identical. For that channel, one does not need to 
employ angular momentum theory, and the Kraus operators are the generators 
of the su(n) Lie algebra [14] . But please be aware that this is not the case for 
our arbitrary spin coupled to a freely fluctuating impurity spin where all higher 
state multipoles K > 1 come with their own different decay rates ([46]) . For 
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this channel, the author has not been able to determine the Kraus operators for 
finite times (but would certainly be happy to receive any valuable information 
on that point by his readers). In this case, there seems to be no easier way to 
describe the spin relaxation than in terms of the irreducible components: 

^) = E e ~ 7if W0)T«f. (52) 

KQ 

This implies the simple exponential decay 

A KQ (t) = e-^ Kt A KQ (0) (53) 
for the irreducible components (j3"5|) of any observable A. 

3 Diffusion 

3.1 Transport 

We may call "transport" a movement from a point r to a point r' that is induced 
by an external cause. In free space, propagation is ballistic: the average square 
of the distance covered after a time t scales like (r 2 ) ~ v 2 t 2 where v is the 
particle's velocity. A disordered medium contains impurities that interrupt the 
ballistic movement. So-called "quenched" disorder is fixed for each realisation 
of an experiment, but varies from experiment to experiment when samples are 
changed. Predictions about observables will involve an average (...) by integrat- 
ing over a classical disorder distribution or by tracing out uncontrolled quantum 
degrees of freedom. Generically, the averaged expectation value behaves as a 
diffusive quantity: (r 2 ) ~ 2Dt with D the diffusion constant. 

In a hydro dynamic description, diffusion is a direct consequence of two very 
basic hypotheses: 

(i) a local conservation law dtn + V • j = 0, also known as the continuity 
equation, linking the local density n(r,t) and the local current density 
j(r,t), and 

(ii) a linear response relation j = -DVn, known as Fourier's law in the 
context of heat transport and Fick's law in the context of particle diffusion. 

By inserting the second relation into the first, one finds immediately the diffusion 
equation (d t — _DV 2 )ra = 0. The hydrodynamic description is only valid for times 
and distances large compared to the scales on which microscopic scattering takes 
place. The linear response coefficient D has to be determined microscopically. 
In essence, the simplest physical process leading to diffusion is a random walk 
or repeated elastic scattering. A kinetic description where point particles collide 
with obstacles permits to derive the diffusion constant D associated with this 
process as function of microscopic scattering parameters. In this section, we will 
derive the appropriate kinetic equation for elastic momentum scattering from 
first principles using a master equation approach and determine the relevant 
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diffusion constant. A largely equivalent presentation with more details can be 
found in chapter 4 of the recommendable book "Quantum Transport Theory" 
by J. Rammer 26] . 

3.2 Momentum - a primer on translations 

If spin is to be understood by considering rotations, then momentum is under- 
stood by considering translations. 

Let a particle be prepared at a point ro = with a spreading described by a 
function f(r). After translation of the entire preparation apparatus by a vector 
a G R d , the new position is described by the function 

[T a f](r):=f(r-a) (54) 

such that the particle is indeed centered around r' = a. In this so-called "ac- 
tive" formulation of translations, relation (|54|) defines the action of the trans- 
lation operator T a in a functional space of, say, probability distributions, in 
exact analogy to the case of rotations treated in section 12.1.31 Here, it is an 
infinite-dimensional representation of the group (R, +) of real numbers with the 
addition "+" as a group law; to be precise, in d dimensions it is the d-fold direct 
product of such representations. Translations have the identity element E = To 
and inverse T" 1 = T_ a . This group is Abelian because different translations 
commute: T a Tb = T a+ f, = Tb+ a = TbT a . Furthermore, this group is a simply 
connected Lie group, and all translations can be generated by exponentiation 
T a = cxp{— ia ■ T} of i = 1, . . . , d generators Ti that form a Lie algebraQ 

In the functional representation (fS~4")) . all translations T a — exp{— ia • T} are 
generated by T = — iV, the derivative with respect to the position coordinate, 
which becomes apparent through a Taylor series expansion 

f(r -a)= f(r) - a ■ V/(r) + • • • = exp{-a ■ V}/(r). (55) 

In quantum physics, the fundamental commutation relation between position 
and momentum, [fi,pj] — iMy, implies that the momentum observable p = hT 
is the translation generator, p — — i?iV in position representation. Translations 
are implemented by the unitary operator 

U{a) = cxp{— ia • p/U} (56) 

such that the particle's position operator r transforms as 

r' = U{a)rU{a)^ = r - a, (57) 

4 If, however, the possible positions lie on a lattice, only discrete translations by lattice 
vectors are allowed, and translations are a representation of (2!,+), the additive group of 
integer numbers. In an infinite volume, the translation group is not compact because it is 
not bounded such that its representation theory is quite different from the compact rotation 
groups SU(n). Notably, it has no finite-dimensional irreducible representations. An exception 
are discrete translations on a lattice with periodic boundary conditions, often used for classifi- 
cation of crystalline lattices. This group is the cyclic group, and its irreducible representations 
are labelled by the admissible wave vectors of the reciprocal lattice |28| . Physicists are taught 
to know this under the epithet of Bloch's theorem. 
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quite similarly to the corresponding identity @ for rotations. In contrast to the 
angular momentum commutation relations ([6]), the simpler version [pi,pj] = 
for the translation generators reflects their Abelian structure. 

The momentum operator is diagonal in the basis of momentum eigenstates, 
p\p) = p\p)- Any operator O is translation invariant if and only if it com- 
mutes with the momentum operator, [0,p] = 0. Therefore, it is diagonal 
in the momentum representation, {p'\0\p) = S pp iOp. Notably, in absence of 
any external perturbation, the Hamiltonian Hq should be translation invariant, 
Ho\p) = £p\p), a property sometimes referred to by the expression "the pt's are 
good quantum numbers" . Henceforth, we will choose units such that K = 1 and 
thus drop the distinction between momentum and wave vectors: p — hk = k. 

The spatial form of the wave functions tp P {r) = (r\p) is fixed by their transla- 
tional properties: using ([54]) together with ([56]) yields (r — a\p) = (r\[U(a)\p)] — 
cxp{— ia ■ p}(r\p) such that (r\p) = C expjir • p] up to a normalisation factor. 
In a finite volume f2 = L d , these plane waves are square- integrable and can be 
properly normalised. Instead, we choose to work in the limit Q — > oo and fix 
C = 1. The identity is then resolved by 

1 = J d d r\r)(r\= J ^\p)(p\- (58) 

Using the plane-wave expansion, Fourier transformation is written using the 
convention 

f{x) = J J^ eiX Pfp With f P = j^xe-^fix). (59) 

This choice is convenient because factors of 2tt are always associated with p- 
integrals which, if required, can be easily converted back to finite- volume sums, 

I^F(p) = L- d E p F(p)- 

3.3 Master equation approach to diffusion 

We will show that a microscopic quantum derivation of diffusive behaviour is 
possible starting from the single-particle Hamiltonian 

H = H a + V. (60) 

Here the translation-invariant part Hq describes free propagation in momen- 
tum eigenstates \p) with eigenenergies e p = p 2 /2m; the generalisation to more 
general dispersion relations is straightforward. The random impurity potential 

V = - x i) = Yj e~ i *"" i »(r)e 1 *-" 1 (61) 

i i 

describes momentum scattering by an arbitrary potential v(r), typically quite 
short-ranged, centered on random classical positions {a^}. Deriving observable 
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quantities will involve averages over all possible realisations of the disorder. The 
ensemble average of any quantity 0({xi]) is the integral 

= J (Hd^Pdx^Oaxi}) (62) 

i 

over all inpurity positions weighted by their distribution P({xi]). The sim- 
plest distribution P({xi}) = Y\, L P{xi) = Yli=i ^ 1 describes N uncorrelated 
impurities with average uniform density n = N/Cl in any finite volume SI. This 
distribution will be used in the following, and we will go to the thermodynamic 
limit N, fi — > oo with fixed density n. 

This description is valid if the mass M of scattering impurities is much larger 
than the mass m of scattered particles. This is the case for electrons scattered 
by lattice defects in solid state devices or for photons scattered by cold atoms 
when recoil can be neglected. The impurities then have no internal dynamics 
and simply realise an external potential v(r); this description can be obtained 
in the limit m/M — > from the more general model, where also the dynamics 
of impurities is taken into account. Note that this limit is just the opposite of 
the usual picture used for quantum Brownian motion (cf. Sec. 3.4 (c) in [TJ), 
where one tracks the movement of a large test particle bombarded frequently 
by smaller ones. 



3.3.1 Derivation of master equation 

Here, we follow the standard derivation of a master equation for open quantum 
systems []j|J [37] by adapting the Born-Markov or weak-coupling recipe to our 
case, cf. Sec. 4.1 in [1!. The ensemble average (1B"2")) now plays the role of a trace 
over bath variables. 

Starting from the Liouville-von Neumann equation dtp{t) = — i[V"(t), 
in the interaction representation A(t) = U~Q(t)AUo(t) and developing to second 
order in the interaction leads to the pre-master equation for the averaged density 
matrix p({xi},t) =: p(t): 

t 



d t p(t) = -i[V(t), P (0)] - / V(t), [V(t - t>), p({xi}, t - t')} dt'. (63) 



We may assume that the initial density matrix p(0) does not depend on the 
disorder configuration - it may for instance represent an initial wave-packet 
prepared far from impurities whose temporal evolution we wish to follow. Then, 
the first term [V(t),p(0)} vanishes because 



V(t) = V(t) = V (64) 

is a constant real number that shifts all energy levels e p of Hq by a uniform 
amount and thus gives no contribution under the commutator. 

The resulting equation is still exact, but not useful: it is not a closed equation 
for the ensemble-averaged p(t) since the density matrix inside the integral still 



21 



depends on the disorder configuration. In order to cope with the integrand, one 
typically proceeds with the so-called Born approximation, replacing the exact 
density matrix inside the integral by its average: 



3 t p{t) = - / V(t), [V(t - t'),P(t - t')} df . (65) 

Now, we are left with an effective Gaussian model of disorder since everything 
depends on the pair correlations VV. Concerning the time dependence, we still 
face a difficult integro-differential equation for pit) . If the time scale of scattering 
is much smaller than the time scale of evolution we are interested in, we can 
perform the Markov approximation by replacing p{t — f) i— > p(t) inside the 
integral and by letting the upper limit t of integration go to oo such that now 
we have a closed differential equation for pit). Reverting to the Schrodinger 
representation we find the following master equation for scattering by fixed 
impurities: 

d t p{t) = -i[H ,p(t)]+Vp(t) (66) 
with a scattering superoperator T> defined by 



Vp{t) = Vp(t)W + Wp(t)V - VWpit) - p(t)WV (67) 



where 



/•oo />0 

W:= Vi-t')dl/ = / ulit')VU (t')dt' =:V]e-**%(f)e*"». (68) 

JO J~oc j 

It will become apparent in section 13.3.51 that the weak coupling or Born ap- 
proximation discards genuine quantum corrections and entails purely classical 
dynamics. Instead of "approximation" , we had better speak of "simplification" 
because at this stage we have no means of knowing whether the resulting de- 
scription is truly an approximation or perhaps qualitatively wrong. And really, 
in section 14.2.11 we will see that quantum corrections need to be considered in 
phase-coherent samples. 



3.3.2 Momentum representation 

In order to see what kind of evolution the master equation (|66[) describes, we 
evaluate it in the momentum representation in which the free Hamiltonian Hq 
is diagonal. In the short-hand notation |1) = |p x ) and E\ — e Pl , the first 
scattering contribution reads 

{l\Vpit)W\A) = ^^(l|e- i P^>w(f)e i P-^|2)(3|e- i P-^u;(f)e i P-^|4)(2|p(i)|3). 

2,3 i.j 

(69) 

2 

The terms with i ^ j are proportional to V and cancel with an equivalent 
contribution in (|67|) that comes with a minus sign (as before, the average V 
gives no contribution thanks to the commutator structure of the equation of 
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motion). In the terms i — j, we can take the translation operators outside the 
matrix elements and perform the ensemble average: 



pj r 

i ^ J 

= n{2ir) d S(p 1 +p 3 -p 2 -p 4 ). 



(70) 



As expected, the average over a uniform distribution restores translational in- 
variance which is equivalent to the conservation of total momentum expressed 
by (2n) d S{p 1 + p 3 - p 2 - p 4 ) =: S 1+3 ^ 2+4 . 

Proceeding in the evaluation of (|69|) , the first matrix element 



(l\v(f)\2) = / d d re- i{p i- p * )r v(r) = « Pi _ Pa =: v 12 (71) 



is the Fourier transform of the scattering potential. In the second matrix element 
(3|w(r)|4) = / oo (3|C/ t (i')w(f)C/ (^)l 4 >di', we can pull out the time integration 

o : 

e i(e 3 -e 4 )t dt > = = . r 34 _ ( 72 ) 

-oo £4 - £3 + lO 

Readers with a background in perturbation theory will recognise this as the 
matrix element r 34 = i(p 3 |Go '(£4)^3) of the free retarded resolvent operator 
G^(w) = (u> — Hq + iQ)^ 1 . This gives us a hint on the applicability of the Markov 
approximation: the rapid time scale here is the inverse energy difference £3 — £4 
of incident and scattered state. The effective evolution on much longer time 
scales into a new state \p 3 ) is constrained by (fT2"|) to the energy shell £4 of the 
incident state \p 4 ). The imaginary contribution of T produces the Lamb shift 
that renormalises the original energy levels, cf. Sec. 4.1 in [T], whereas the real 
part yields the relaxation rates that render the dynamics irreversible. 
Altogether, this first contribution to the collision functional reads 

(l\Vp(t)W\4) = nJ2 S i+3,2+i^v 3 4T 34 {2\p(t)\3}. (73) 
2,3 

Collecting all four contributions gives 



(l|X>p(i)|4) - nJ2 [<5i+3,2+4«i2f34(ri2 +r 34 )(2|p(i)|3) 

2,3 

-<MM 2 r 2 i + M 2 r 43 )(i|,5(i)|4) 

We choose to use the parametrisation 

Pi=P + <7/2, p 2 =p' + q/2 
Pi=P 9/2, p 3 = p' - q/2 



(74) 



(75) 
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that complies with the conservation of total momentum ([70]) . With this parametri- 
sation, V\2 — v p ~ P ' — v^ 4 . Furthermore, energy differences become 

(p-p')-q , . 

El - e 2 =: £ P - s P > H s ( 76 ) 

2m 

and £1 — 64 = p - q/m, and we write the sum of matrix elements (|72[) in the form 

ri 2 + r 34 =:rv(g). (77) 

Finally, disposing of the overbar p(t) t— > p(t), the density matrix elements 
p(p, q, t) :— (p + %\p(t)\p — % ) obey the master equation 

d t p(p, q, t) = -i— q, t) + V[p{p, q, t)] (78) 



t)}. (79) 



with the scattering functional 

/d d p' 
j2^d\ v p-p'\ 2 { t pp'(<i)p(p' - r P 'p(q)p(p,q, 

3.3.3 Trace conservation and continuity equation 

Clearly, the master equation ([75)) has the form of a kinetic balance equation 
where the scattering functional (|T9"|) contains transitions p' — > p that increase 
the magnitude of p(p, q, £), and also negative contributions of depleting transi- 
tions p p' . Since we have neither sinks nor external sources, the net effect 
must be zero which should be apparent in a conservation of the local proba- 
bility density. And really, the master equation first of all preserves the trace, 
d t tr{p(t)} — J ^Mp d t p{p, 0, t) — 0, since the collision functional is antisymmet- 
ric under the exchange p <-> p' and thus 

/ (^ P W^)]=0. (80) 

This antisymmetry is inherited from the double-commutator structure of (|65p 
and holds for all spatial Fourier momenta q. Summing the q-dependent master 
equation over p thus leads to the continuity equation 

d t n q (t)+iq-j q {t) = (81) 

that links the Fourier transforms of the local density n(r,t) and local current 
density j(r,t), given as the first two p-moments of the density distribution: 

n q {t) = J d d re-^ r n(r,t) = J APp(p, q ,t), (82a) 

j q (t) = J d d re-^j{r,t) = J ^E-p(p,q,t). (82b) 

The current vanishes by parity for isotropic distributions p(p, q, t). 
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3.3.4 Momentum isotropisation 

What kind of dynamics does the master equation describe? A first, simple an- 
swer is possible by considering the limit q = that describes spatially averaged 
quantities. The definitions (|T2"|) , (|7rJ|) and ([77)) imply r pp '(0) = 2n5(e p — e p >) 
which assures the conservation of energy during elastic scattering. Since the 
isotropic energy e p fixes the modulus of p' = ph , the angular probability dis- 
tribution p(pn,Q,t) =: f e (n,t) at fixed energy satisfies 

d t f e (n,t) = 2nn J ^S(s - e P 0K- P '| 2 [fe(n',t) - f s {n, t)] . (83) 

For an isotropic point-scatterer potential v(r) = vo5(r) that has no dependence 
on momentum, the equation of motion takes the simple form 

dtf e (n,t) = - 7 ci(e) [fe(n,t) - (f e (t))] (84) 

where (f 6 (t)) ■= J dil' d f e (n',t) is the angular average of the distribution, prop- 
erly normalised: d£l' d is dO' /2tt in d = 2 dimensions and d^>'d(cos6 |/ )/47r in 
d = 3. The elastic scattering rate 

7ci(e) - 2w{e)nvl (85) 

is defined in terms of the density of states v{e) = J ^ 2 ^ d 5(e—e p ') and can equally 
well be obtained by Fermi's golden rule. Clearly, equation describes a simple 
exponential decay of the initial angular distribution f £ (n 1 t = 0) towards a 
completely isotropic distribution (f £ (t)). Thus, our model of elastic momentum 
scattering by fixed impurities leaves the kinetic energy conserved, but describes 
the isotropisation of the momentum distribution with a rate 7 e i. 

By the same token, a global net current j (t) initially different from zero 
decreases to zero exponentially fast. In other words, an initial wave packet 
launched with a definite velocity loses the memory of its initial direction on a 
time scale t q \ — l/7 i- 

3.3.5 Boltzmann-Lorentz equation 

In order to know on what spatial scale the momentum isotropisation occurs and 
how the average position of the wave packet evolves in time, we have to consider 
the master equation ((78|) at finite q. We expect diffusive behaviour to appear on 
large spatial scales and thus make a Taylor expansion around q = 0. Retaining 
only lowest-order terms in q/p <C 1 (which corresponds to spatial scales much 
larger than the particle's wave length), the equation of motion becomes 

d t p(p, q, t) + i — p(p, g, t) = C[p(p, q, t)} (86) 
m 

with the purely elastic collision integral 
f d d n' 

C[p{p, q, t)] - 27m / S(e p - e p ,)\v p - p > | 2 [p{p', q, t) - p(p, q, t)] . (87) 
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The only explicit occurrence of q in ([56"]) is now in the ballistic term on the 
left-hand side that originates from the free evolution with Hg. 

The parametrisation (|75|) has the additional advantage that the density ma- 
trix elements p(p, q, t) are the spatial Fourier transform of the Wigner distribu- 
tion [25] 

W[P ' r ' t] = JTW I & d r'{r- r -\p{t)\r+ r -yv- r '/ n 

a (88) 
1 f d d q ■ 

that represents the quantum density operator in a classical phase space. Here, 
we have momentarily restored fi.'s visibility. With this standard normalisation, 
the probability density n(r) = (r\p\r) with normalisation J d d rn(r) = 1 is given 
as the marginal n(r) — J d d pW(p,r). Conversely, the momentum distribution 
is w(p) = J d d rW(p,r) = (2nh)- d p(jp,0), and / d d rd d pW{p,r) = 1. 

A Fourier transform with respect to q now yields the kinetic equation 

d t W(p, r, t) + —p- V r W(p, r, t) = C[W(p, r, t)} (89) 
m 

with the same elastic scattering integral (|87l) . Remarkably, we have obtained 
precisely the linear Boltzmann-Lorentz equation for the classical phase-space 
density W(p,r,t) under elastic scattering from fixed impurities [3D]. For pho- 
tons, this equation is known as the radiative transfer equation. Ex post, we can 
therefore conclude that the Born approximation ([65]) was the crucial step that 
discarded quantum corrections to propagation amplitudes and left us with a clas- 
sical phase-space distribution. This interpretation transpires also by analysing 
Feynman diagrams in perturbation theory and path-integral approaches (see 
chapter 4 of [5S] for details). 

The weak disorder limit in which the Born approximation and Boltzmann 
transport theory are expected to be valid corresponds to the regime where dis- 
order corrections to the free energy of the particle are small, 7 e i(e) *C s. The 
scattering time defines a typical length scale, the elastic scattering mean free 
path £ e i = vqt c \ or average distance between successive scattering events. The 
weak disorder limit can also be stated as l/{ki c \) <C 1 which requires that suc- 
cessive scatterers are placed in the scattering far field: the average distance £ e \ 
must be larger than the wavelength A = 2n/k. This is a low-density argument 
because £ e \ — l/(na c \) in terms of the density n of scatterers and their total 
elastic scattering cross section a c \. 

It is allowed to neglect the explicit ^-dependence inside the collision integral 
(|87[) . which contributes already a factor j e i, on hydrodynamic scales q£ c \ <C 1. 
Consistently, it is precisely in the hydrodynamic regime that we wish to deter- 
mine the diffusion constant. We have already derived the continuity equation 
(j8T]) . Making step (ii) of the general argument presented in 13. 1( we now turn to 
the calculation of the diffusion coefficient in the linear response regime. 
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3.4 Linear respose and diffusion constant 



Any phase-space distribution W cq (p) that is homogeneous, independent of time 
and rotation invariant, i.e., depends only on the modulus of p, is a solution of the 
Boltzmann-Lorentz equation (|89[) since each term vanishes separately. The cor- 
responding density matrix elements are of the form p eq (p, q, t) — (2ir) d 5(q)p cq (p). 
This type of solution is called a global equilibrium. The underlying statistics 
could be a Fermi-Dirac or Bose distribution, or their classical limit, the Boltz- 
mann distribution. 

By creating a small gradient of concentration, one can then induce a linear- 
response current that is proportional to the driving gradient; the coefficient 
of proportionality is the diffusion constant. Kinetic theory permits to calcu- 
late linear response coefficients. We will follow the linearisation method devel- 
oped by Chapman and Enskog [3D] m order to derive the diffusion constant, 
but in terms of the density matrix components p(p, q, t) instead of the space- 
dependent Wigner distribution because it will prove useful to work with Fourier- 
transformed quantities. 

3.4.1 Linearisation a. la Chapman-Enskog 

Suppose that initially, the distribution function p{p 1 q, 0) is a local equilibrium 
solution (i.e., isotropic in momentum) established by local scattering on a rapid 
time scale r cl = However, we assume a non-delta like dependence on q, 

i.e., a finite gradient in real space. It is therefore no longer a global equilibrium 
solution of the Boltzmann-Lorentz equation. The time r cq it takes to reach 
global equilibrium is much longer than the local scattering time such that we 
have a small parameter T e i/T eq <C 1 (sometimes called the "Knudsen number"). 
The linearisation method of Chapman and Enskog works by expanding the 
distribution function formally in powers of this small parameter, 



For the linear response calculation, these first two terms will suffice. The col- 
lision integral effectively multiplies the distribution by 7 c i = 1/t c i such that 



which is satisfied if po(p, q, t) is locally isotropic. By parity, the current density 
([82b]) then is entirely generated by the correction, 



The continuity equation (15T]) then implies that to lowest order the local 
density remains time-independent: dtPo(p,r,t) = 0. To this order = 1, the 
kinetic equation reduces to 



p{p,Q,t) = Po(p,q,t) +r e ipi(p,r,t) + 0(r e 2 i). 



(90) 





(92) 



q p 



Po(p,q) =C[T clPl (p,q)}. 



(93) 



m 
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In order to calculate the current (|9"2"|) , we need to solve this equation for pi as 
function of po- We content ourselves with the simple case of an isotropic point- 
scattering potential v p — v treated in section 13.3.41 (other cases can be treated 
by an expansion in angular eigenfunctions, see [3T] for the case of a potential 
with scattering anisotropy). Then, using the right-hand side of (|54"|) , we find 

C[T el p 1 ] = -p 1 (p,q) (94) 

plus an isotropic term (pi) that does not contribute to the current (|92|) anyway 
and can be dropped such that 

pi{p,q) = -i — po(p,q)- (95) 



3.4.2 Diffusion coefficient 

Inserting (|9"5"|) into (|9"2"|) . the resulting current reads 

i f d d p 

U = -^2 J ^pP(9-P) T ci( e p)Po(p,g). (96) 

By isotropy, J d d ppiPjf(p) — (Sij/d) J d d pp 2 f(p), and the current is collinear 
with q. To lowest order in q we find 

j q = -iqD n q , (97) 

the Fourier-transformed version of the linear response relation j = —DoS7n. 
Kinetic theory has allowed us to to calculate the diffusion coefficient 

D = ^ (98) 

as the product of an effective velocity and scattering time averaged over the 
momentum distribution, 

1 f d d p 

^= ^2 J j^dP 2 T c i(e p )p {p). (99) 

Often, the distribution po(p) is a sharply peaked function around a certain 
momentum po (for instance the Fermi momentum pf for electrons), whereas 
p 2 r e i(e p ) according to ([83)1 varies smoothly with the density of states such that 
vfjr = pQT o \(e P0 )/m 2 . In terms of the scattering mean-free path £ e \, the diffusion 
constant for isotropic point scatteres can also be written 

Do = (100) 

For anisotropic scattering, one has to replace the scattering mean-free path £ e \ 
by the transport mean- free path £ tT [3Tj . 
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3.4.3 Diffusion 



Inserting the linear response current (|9T|) in the continuity equation leaves us 
with a simple differential equation d t n q {t) = —DQq 2 n q {i) that is immediately 
integrated to give an exponential decay 

n q (t) = c- D °« 2t n q (0) (101) 

for the Fourier components of the initial density fluctuation. Long-range fluc- 
tuations (q — > 0) take a very long characteristic time r q — l/(Doq 2 ) — > oo to 
relax because of the constraint imposed by local conservation. This diagonal 
decomposition into Fourier modes with their continuous momentum index q is 
the analog of the spin relaxation (|55|) where the discrete index K separates high- 
er irreducible modes with rather large relaxation rates 7^ from the isotropic 
component K = or trace that is conserved. 

The expectation value of the average radius squared, 

= tr{r 2 p(t)} = -V^(t)| q=0) (102) 

for the diffuse density (|101|) reads: 

F)(*) = (r 2 ) + 2<tD t. ( 103 ) 

As expected, the long-time behavior of the particle's displacement is indeed 
governed by the Boltzmann diffusion constant (|100|) . 

4 Diffusive spin transport 

Having described in section [2] the relaxation of a single motionless spin in a 
fluctuating field, and in section [3] the diffusion of spin- less massive particles, we 
now combine these two pictures and consider a spin on a massive carrier particle 
that moves and encounters impurities. The Hamiltonian is still of the form 
H = Hq + V where Hq describes ballistic propagation in momentum and spin 
eigenstates \pa) :— \p) <8> \sa) with spin-independent eigenenergies e p — p 2 /2m 
(spin quantum numbers will from now be called a in order to avoid confusion 
with the particle's mass). The impurity potential V could describe momentum 
scattering, spin-flip scattering, spin-orbit coupling, and the like, by randomly 
distributed scatterers. 

In the following, we will consider in detail the case of freely orientable mag- 
netic impurities that induce spin-flips. Other mechanisms such as spin-orbit 
scattering can be treated along the same lines. Actual laboratory realisations 
include electronic spin-flip scattering, quite relevant even for very low impurity 
concentrations, and the randomisation of photon polarisation under the influ- 
ence of scattering by atoms with degenerate dipole transitions in cold atomic 
clouds. 

5 This analogy becomes even clearer if positions in a finite volume are restricted to a lattice 
of L d sites, because the irreducible representations of the discrete translation group (which 
under periodic boundary conditions is the d-fold direct product of the cyclic group 'El) are 
precisely labelled by the different allowed q- vectors of the reciprocal lattice |28| . 
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4.1 Master equation approach to diffusive spin transport 
4.1.1 Deriving the master equation 

In addition to the elastic momentum scattering potential (f6Tj) that for clarity 
we now call V e i, consider then a spin- flip interaction potential 

V tt = Y^v tt (r-x j )S-T j (104) 

3 

between the spin S and a collection of freely orientable magnetic impurities 
modelled as spin ^ with Pauli matrices Tj centered at sites Xj. The mag- 
nitude of spin-spin coupling has been included in the short-ranged potential 
v s f(r) whose spatial dependence induces momentum scattering. The ensemble 
average now contains the usual average (|62|) over random sites as well as a 
trace ti{ T y(p{ T y) over impurity spins. We will assume that these impurities are 
distributed independently and isotropically, p{ T j = • pj with pj = ^Ig. 

It is now a simple task to derive a spin-diffusion master equation for the 
density matrix p(t) that operates on the combined Hilbert space TL = Ti s ® TL P 
of spin and momentum by retracing exactly the steps of sections 13.3.11 and 
13.3.21 with the new potential V — V c \ + V a { instead of just V c \. Within the 
Born approximation second order in V, mixed terms T4fT4i = V s { V e i give no 
contribution (just as the product of averages before), and we can consider the 
impact of I^f separately. A typical term arising in the new spin-flip part is the 
counterpart of (|69|) . (l|V£fp(t)Wsf|4). Together with the momentum-scattering 
factors appearing already in (|T3[) . we find now an additional sum over spin 
indices that defines the action of a spin-flip superoperator 

(oi|Vp|o- 4 ) = ^2 {°i\ s ■ T^ie^S ■ t\cj a ){(J2\p\v?,) 

(105) 



where the isotropic average over the impurity spin leads to the scalar contraction 
itr T {rjT 7 } = Sij (cf. exercise 1). 

Collecting all terms then gives the master equation that describes elastic and 
spin- flip scattering in the small-q limit: 

d t p(p,q,t)+i^p(p,q,t)=C[p]+£[p] (106) 

Here, the elastic collision integral 

C[p] = ^rj ~ V) W, ^ t) - ^ *)] (107) 

describes momentum isotropisation with a total scattering rate 

7 = 7el + Isi (108) 
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that includes the spin-flip contribution 7 s f = s(s + l)2TTi/(e)v^ i n s {. Summing 
the rates, a prescription known as "Mathiessen's rule" [32], is permitted if the 
scattering mechanisms do not interfere with each other which is the case in the 
low-density Born approximation considered here. Note that the collision integral 
(|107|) acts solely on the momentum degrees of freedom, but is the identity in 
spin space. Genuine spin-flips are generated by 

£\P\ = 4t / WTi S ^~ e P ')Cp(p',q,t) (109) 
v(e) J {2ir) d 

where £ — "f s f£ is precisely the spin relaxation Lindbladian (|24[) derived in 
section \2. 2. 51 now with the spin- flip relaxation rate 7 s f. 



4.1.2 Diffusive spin relaxation 

Solving the spin-flip master equation (|106[) in the linear response regime is now 
exceedingly simple by projecting it onto irreducible spin components p( K > — 
PkqTq since the spin relaxation Lindbladian is diagonal in that basis, 
£p( K ) = \ K p^ K \ with eigenvalues Xk given by (|4l)j). 

Summing the master equation for the Kth spin sector over p gives the con- 
tinuity equation 

ftnW(t) + iq ■ = -iKn^Ht) (110) 

for the density = J p( K ) (p, q) and associated current density ■ 

The spin-flip Lindbladian is responsible for the appearance of a source term on 
the right hand side, or rather a sink with spin relaxation rate 

K(K+l) 

7A- = |A A -|=7sf 2a{a + 1 y (HI) 

In the limit q — ► 0, we recover exactly the global spin relaxation (|48p of section 
12.3.61 In particular, the total trace is conserved since the spin trace is the K = 
sector with vanishing eigenvalue Ao = 0. 

In the linear response regime, the Chapman-Enskog method of section [3.4.1l 
carries through in each spin sector K . A difference occurs for the time deriva- 
tive of the locally isotropic component py K \p,q,t). The master equation (|106p 
implies 

d t pi K) (p,q,t) = ~-tKP ( K \p,q,t) (112) 

which is solved by p$ (p, q, t) = c tKt po(p, q, 0) instead of just being con- 
stant as in 13.4. 11 Finally, the spin-sector density components show diffusive as 
well as spin-relaxation dynamics 

nf\t) = e- D °« i '*e-' ), **nW(0) (113) 

where the Boltzmann diffusion constant Dq is evaluated with the total momen- 
tum relaxation time r = 1/7 from (|108p . 
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Now we can answer the question that was the starting point of the lecture 
(cf. Fig. [I]): Imagine that we can inject spin-polarised particles |f) := \s, a = +s) 
with probability p| (0) = (T|po|T) = 1 on one end of a diffusive medium of length 
L. What is the probability pj (L) of retaining the spin polarisation at the other 
end, assuming that we have spin-sensitive detection? 

By taking the matrix elements (T|/° (*)|T) of each irreducible spin com- 
ponent, we find that the probability relaxes during the transmission time t = 
L 2 /2D a as 

1 "3q 2 

Mj 2.s + l + S (s + l)(2 S + l) e +•■■ (U4) 
= i(l + e- t ^) 

with the last line valid for electrons for which 1/t\ = 47 s f/3. Equivalently, the 
degree of spin polarisation 

«*) = M %Tm ='- t/T1 (115) 

simply relaxes from unity to zero. Naturally, for long enough times, the disti- 
bution relaxes towards its equilibrium value p cq = 2s 1 +i t° nave an y magnetic 
quantum number a = — s, . . . , s. 

In terms of length scales, the spin relaxation time permits to define the spin 
relaxation length A s f = ^/2D r s f. In solid-state devices, the density of magnetic 
impurities can be controlled such that 7 s f <C 7el which means that spin coherence 
can be maintained quite efficiently, even on scales L 3> t e \ where the momentum 
dynamics is no longer ballistic, but already diffusive. 



4.2 Quantum corrections 

In the Boltzmann transport theory developed in section [31 one propagates ef- 
fectively classical probabilities. However, in an environment that preserves the 
phase coherence of the wave, one must propagate probability amplitudes that, 
by the superposition principle, allow for interference phenomena. Elastic impu- 
rity scattering does preserve the phase coherence of the propagating wave which 
means that we have to expect quantum corrections to the Boltzmann transport 
theory whenever external phase-breaking mechanism are so rare that the cor- 
responding dephasing time is much longer than the elastic scattering time 

Tel. 



4.2.1 Weak localisation 

A prominent example for such a quantum correction is weak localisation (WL): 
the resistance of weakly disordered metallic samples shows a negative magne- 
toresistance dp/dB > at small fields if spin-orbit scattering is absent [33l [34] , 
This contradicts the classical Boltzmann-Drude picture that predicts that the 
resistance should increase with a magnetic field. The reason can be understood 
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Figure 3: Diffusive path (above) with the classical Boltzmann contribution of co- 
propagating amplitudes (left) and the corresponding counter-propagating amplitudes 
(right) that lead to weak localisation corrections. 



by considering the quantum return probability to a point which includes the con- 
structive interference of counter-propagating amplitudes and is therefore larger 
than the classical probability, see figure [31 The interference effect is masked 
by a large enough magnetic field since loops of different sizes pick up different 
Aharonov-Bohm phases. 

Without any dephasing mechanism nor an external magnetic field, the inter- 
ference correction D = Do + AD to the classical diffusion constant Do is given 
by summing the contributions of all closed diffusive paths, 

AD _ 1 f f°° j -D q 2 t [ d d q 1 

^--^M7)JJ^yJo dte ° a i WV' (116) 

In low dimensions d = 1,2, this q- integral diverges for small q which indicates 
that WL corrections become very important whenever the interfering amplitudes 
can explore large length scales. The formal divergence is cured by a cutoff that 
can be either the system size itself or a finite phase-coherence length — 
WDoT<f, due to dephasing mechanisms on a time scale . [£] 

4.2.2 Dephasing of weak localisation by spin-flip scattering 

Including spin-flip scattering into the weak-localisation picture can be done by 
diagonalising the spin- flip vertex appropriately [53]. The irreducible subspaces 
turn out to be the usual singlet and triplet state subspaces spanned by the 
Hilbert-space vectors \KQ)E 

To sketch the result, the WL correction is written 

^ = _^_/^y m (X17) 



°The ultraviolet divergence for large q is cut off by the shortest scale of scattering, t e \, and 
irrelevant for the present considerations. 

7 Not operators \KQ) in Liouville space as for the propagated intensity: the difference is due 
to the fact that the two amplitudes in weak localisation loops propagate in opposite directions 
such that two states have to be recoupled instead of a state and its complex conjugate. 
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where each spin channel comes with a weight wk = (— ) 2s+K (2K + 1)/ (2s + 1). 
For electrons, wo — — h and Wi = |. More importantly, each spin channel is 
damped with its coherence time t c {K). They are given in terms of the spin 
relaxation rates PS)) by a recoupling relation that reduces to 

t c (K) T Bl t s1 V 4s(s + 1) J 

such that for electrons r c (0) = r s f/2 and r c (l) = 3r s f/2. 

Clearly, the interference in all spin channels is dephased rather efficiently 
with a rate given essentially by the spin-flip rate, the numerical prefactor being 
entirely fixed by geometry. This effect is visible even in quite pure samples as 
admirably shown by F. Pierre and coworkers |35j . see figure 2J The complete 
theory needed for linking the experimental data to impurity concentrations re- 
quires to take into account the interaction between electronic excitations near 
the magnetic impurity (Kondo physics) as well as the temperature-dependent 
dynamics of impurity spins (Korringa physics) that have both been neglected 
within the present lecture; for details, see [55] . 

Multiple coherent scattering of photons in clouds of laser-cooled atoms is 
also subject to spin-flip physics: the so-called coherent backscattering effect (an 
interference enhancement of backscattered intensity) is severely reduced when 
photons are scattered by dipole transitions with a Zeeman degeneracy |36j . 
The analytical theory for multiple coherent scattering of polarised photons by 
degenerate dipole transitions employs the concepts of irreducible decompositions 
exposed in section [2j Compared to electrons, the theory appears much simpler 
because at low light intensity, photons do not interact. However, the treatment 
of photon propagation is more involved because the field transversality adds 
another source of polarisation relaxation that needs to be taken into account 

Thus, magnetic impurities are a very efficient source of dephasing for interfer- 
ence of spin-carrying particles: the large ground-state degeneracy implied by the 
random orientations of freely orientable impurity spins permits dephasing even 
at zero temperature - when other decoherence processes like electron-phonon or 
electron-electron scattering are suppressed - because no energy exchange is in- 
volved and stocking which-path information comes for free. In return, whenever 
the impurity degrees of freedom can be constrained by other means, then per- 
fect coherence is restored. This has been shown in Aharonov-Bohm interference 
experiments with electronic samples subject to a strong external magnetic field 
that aligns the impurity spins [38, 39 . Similarly, in atomic clouds an external 
magnetic fields that lifts the internal atomic Zeeman degeneracy can be used to 
enhance the effective phase coherence length of diffusing photons [40] , 
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FIG. 5. Phase coherence time as function of temperature in sev- 
eral silver wires. Sample Ag(6N)c (•) is made of the purest silver 
source. Samples Ag(5N)b (O), Ag(5N)c Mn03 (□), and 
Ag(5N)d Mnl ( ) were evaporated simultaneously using our 5N 
silver source. Afterward, 0.3 ppm and 1 ppm of manganese was 
added by ion implantation respectively in samples Ag(5N)c Mn03 
and Ag(5N)d Mnl . The presence of very dilute manganese atoms, a 
magnetic impurity of Kondo temperature T K = 40 mK, reduces 
leading to an apparent "saturation" at low temperature. Continuous 
lines are fits of t^(T) taking into account the contributions of 
electron-electron and electron-phonon interactions (dashed line) 
and spin-flip collisions using the concentration c mag of magnetic 
impurity as a fit parameter (dotted line is r sf for c mag = 1 ppm). Best 
fits are obtained using c mag =0.13, 0.39, and 0.96 ppm, respectively, 
for samples Ag(5N)b, Ag(5N)c Mn03 , and Ag(5N)d Mnl , in close 
agreement with the concentrations implanted and consistent with 
the source material purity used. 

Figure 4: This "set of experiments suggests that the frequently observed 'saturation' 
of in weakly disordered metallic thin films can be attributed to spin-flip scattering 
from extremely dilute magnetic impurities, at a level undetectable by other means"; 
reprinted with permission from F. Pierre et al, Phys. Rev. B 68, 0854-13 (2003). 
Copyright (2003) by the American Physical Society. 
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